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Abstract 

For submicron particles shaped as any axisymmetric body and made with standard canted 
antiferromagnet like hematite or iron borate, the ground state may comprise of a magnetic vortex 
with topologically non-trivial distribution of the sublattice magnetization I and planar coreless 
vortex- like structure for the net magnetization M. For antiferromagnetic particles in the vortex 
state, in addition to low-frequency modes, there are high frequency modes with frequencies over 
the range of hundreds gigahertz, including a mode localized in the region of the radius 30-40 nm 
near a vortex core. 

PACS numbers: 75.75.+a, 76.50.+g, 75.30.Ds, 75.10.Hk 
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The magnetic properties of submicron ferromagnetic (FM) particles of the shape of cir- 
cular cylinder (magnetic dots) are today under great consideration mainly due to their po- 
tential applications Circular dots possess the equilibrium magnetic configuration which 
corresponds to a vortex structure just above the single domain length scale, with the ra- 
dius R > Ran ~ 100-200 nm. The FM vortex state consists of an in-plane flux-closure 
magnetization distribution and a central core with radius of order of 20-30 nm magnetized 
perpendicularly to the dot plane. Reducing of the magnetostatic energy comes at the cost 
of a large exchange energy near the vortex core, as well as the magnetostatic energy caused 
by the core magnetization. For data storage purposes, vortex state is of considerable in- 
terest, because the low magnetic stray field reduce the particles interaction and leads to a 
high magnetic stability of the data written. The magnon modes for dots in virtue of effects 
of spatial quantization possess a discrete spectrum. The possibility to manage localization 
and interference of magnons spawned an idea of so-called magnonics based on usage of these 
modes for developing a new generation of microwave devices with submicron active elements 

n 

|2J. These particles also provide an ideal experimental system for studying static and espe- 
cially dynamic properties of relatively simple topologically non-trivial magnetic structures 
which are fundamentally interesting objects in the research area of magnetism. 

All previous studies of magnetic vortices caused by magnetic dipole interaction were car- 
ried out on the magnetic particles made with soft FM with high magnetization M s like 
permalloy with 4ttM s ~ 1 T. In this Letter we have shown that the vortices can be a 
ground state for sub-micron particles made with another important class of magnetic mate- 
rials, antiferromagnets (AFM), with easy-plane (EP) anisotropy and Dzyaloshinskii-Morya 
interaction (DMI). 

For AFM, exchange interaction between neighboring spins facilitates antiparallel spin ori- 
entation, which leads to the structure with two antiparallel magnetic sub-lattices, Mi and 
M2, \Mi\ = \M%\ = Mq. As typical AFM we can mention hematite a-Fe203, iron borate 
FeBC>3, and orthoferrites, see [3| . These materials are characterized by high temperatures 
of magnetic ordering and have unique physical properties: orthoferrites and iron borate are 
transparent in optical range and have a strong Faraday effect, the magnetoelastic coupling 
is quite high in hematite and iron borate . They possess small but non-zero net magneti- 
zation caused by a weak non-collinearity of sublattices (sublattice canting) originated from 
DMI. 
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To describe the structure of AFM, it is convenient to introduce irreducible combinations 
of the vectors M\ and M, the net magnetization M = Mi + M2 = 1M§m and the vector of 
sublattice magnetization, / = (Mi — M 2 )/2M . The vectors m and I are subject to constraint 
{m-l) = 0, m 2 + P = 1. As \m\ <C 1, the vector / could be considered as a unit vector. The 
mutual orientation of sublattices is determined by a sum of the energy of uniform exchange 
W ey , = H ex M m 2 , and the DMI energy, Wbm = 2M H^(d ■ (m x /)), the unit vector d is 
directed along the symmetry axis of a magnet. The parameters if ex ~ 3 ■ 10 2 — 10 3 T and 
Hy> ~ 10 T are exchange field and DMI field, respectively. Using this energy and dynamical 
equations for Mi and M2, one can find, 

M = M DM (d x i) + — ^ x - j , M DM = ^— , (1) 

where the first term gives the static value of AFM net magnetization Mom, comprising a 
small parameter, Hd/H cx ~ 10~ 2 , second term describes the dynamic canting of sublattices, 
3). Mdm is much smaller than Mq or the value of M s for typical FM, but 



see for details 



the role of the magnetostatic energy caused by Mdm could be essential, and could lead to 
the appearance of a domain structure for AFMs [3j, |4j- We will show that for formation of 
equilibrium vortices in AFM even have some advantage compared with soft FM. 

The dynamical properties of AFM are essentially different comparing with FM. A spin 
dynamics of an AFM can be easily described in the framework of so-called sigma-model 
equation (cr-ME), a dynamical equation for the vector 1 only, with the magnetization M being 
a slave variable [3|. In contrast to the Landau - Lifshitz equation for a FM magnetization, 
the cr-ME contains a dynamical term with a second time derivative of I, combined with 
gradients of I in the Lorentz-invariant form d 2 l/dt 2 — c 2 V 2 Z. For this reason, for AFM 



two magnon branches (instead of one, for FM) exist. The chosen speed c = 7^ AH ex /Mo 
plays roles of both magnon speed and the limit speed of domain walls, it is determined 
by exchange interaction only and attains tens km/s, c ~ 1.4 • 10 4 m/s for iron borate and 
c ~ 2-10 4 m/s for orthoferrites [3]. For both modes the elliptic polarization of the oscillations 



For EP 



of Mi and M 2 is such that the oscillations of the vector I have linear polarization 
AFM these two branches are low-frequency quasi-ferromagnetic (QFM) branch and a high- 
frequency quasi-antiferromagnetic (QAF) branch. QFM magnons involves the oscillations 
of the vectors I and M in the EP, with weak deviation of M from the EP caused by last 
summand in ([1]). The second QAF branch corresponds to the out-of-plane oscillations of I 



with the dispersion law cjqaf(^) = yWg + c 2 k 2 , where k is the magnon wave vector. The 



gap of QAF branch, u g = ^^/2H CX H & contains large value H ex and attains hundreds GHz. 
Thus both magnon frequency and domain wall speed for AFM dynamics, comparing with 
FM, contain a large parameter ^JH CX / H & ~ 30 — 100, H cx and H & are exchange field and 
anisotropy field, respectively, which can be referred as exchange amplification of dynamical 
parameters of AFM. The frequencies of AFM magnon modes u g reaches hundreds GHz, 
with values 170 GHz for hematite, 100-500 GHz for different orthoferrites and 310 GHz for 
iron borate [5J. Recent studies showed a possibility to excite spin oscillations of non-small 
amplitude for orthoferrites [6] and iron borate [7J with the use of ultra-short laser pulses. 

Spin distribution for AFM can be described by the energy functional of the form W[l] + 
W m . Here W[l] describes the energy of non-uniform exchange and the anisotropy energy 
through only the vector I] W m is magnetic dipole energy, 

W m = -\\ MH m al 3 x, (2) 
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where H m is demagnetization field caused by AFM magnetization (TJ|). The sources of H m 
can be considered as formal "magnetic charges" , both volume charges equal to divM and 
surface charges equal — M ■ n, see monographs [J, 18] for general consideration and [91] for 
application to vortices. 

For a pure uniaxial model of an AFM which is applicable for hematite and iron borate, 
and to some extend for orthopherites, W[l] can be presented as, 

W[l] = \J [A{Vlf + K ■ l 2 z ]d 3 x , (3) 

where A is the non-uniform exchange constant and K is anisotropy constant, the xy— plane 
is the EP for spins. Variation of the energy W[l] gives a general two-dimensional (2D) vortex 
solution for the vector I of the form 

— * 

I = e z cos6 + sm9[e x cos(x + y? ) + e y sin(x + ipo)] (4) 
where 9 = 9{r), r and x are polar coordinates in an EP of a magnet, the vector e z is 



the hard axis, the value of tpo is arbitrary, see lOj, lUj. The function 9{r) exponentially 
tends to 7r/2 at r ^> /o, with the characteristic size Iq = a/ 'A/K, and in the center of the 
vortex (at r = 0) sin6 l (0) = 0. Near the vortex core / deviates from the EP which leads 
to the loss of the anisotropy energy. The state (jlj) is non-uniform, what corresponds to the 



loss in the exchange energy. Therefore, for the EP model without taken into account the 
magnetic dipole interaction the appearance of a vortex costs some energy, i.e. the vortex 
corresponds to excited states of AFM. Vortex excitations are important for description of 
thermodynamics of 2D AFM, see lo| . 

For small particles made with canted AFM the energy loss caused by a vortex can be 
compensated by the energy of magnetic dipole interaction. To explain this, note that for a 
uniform distribution state the contribution W m unavoidably results in a loss of the system 
energy, which is proportional to the particle volume V {4, 8]. The energy of the uniform 
state could be estimated as 

^(homog) = 27T NM 2 DM V = 27cNM$(2H D /H e ) 2 V, (5) 
where N is the effective demagnetizing factor in the direction perpendicular to the particle 



axis 4, 8j. In contrast, for the vortex state (jlj) with a chosen value of sin<£>o = 0, with 
M oc (d x T), one can find 

M = a M DM • sin 9(—e x sin x + e y cos x), (6) 

where a = cos<£> = ±1- A unique property of the state (jSj) is that it can also exactly 
minimize the energy of the magnetic dipole interaction W m , giving H m = in the overall 
space. Indeed, the projection of M on the lateral surface of any axisymmetric body with 
the symmetry axis parallel to z-axis, as well as divM, equal to zero. Moreover, in virtue 
of symmetry of DMI (M • d) — 0, the distribution of the magnetization M (J6j) is purely 
planar (in contrast to I) and the out of plane component of M is absent. In the vicinity of 
the vortex core the length of vector M decreases, turning to zero in the vortex center, see. 
Fig. 1. Such feature is well known for domain walls in some orthoferrites jl2 |. Thus the 
AFM vortex is the unique spin configurations which do not create demagnetization field in 
a singly connected body (for FM with \M\ = const a configuration with H m = is possible 
only for a magnetic rings having the topology of torus). 

Let us compare the energies of the vortex state and the uniform state for the AFM particle 
shaped as a cylinder with the height L and the radius R. For the vortex state H m = overall 
the volume of the particle, and the vortex energy is determined by the simple formula 

E y = nAL\n(pr), (7) 



FIG. 1: (Colour online) The structure of the AFM vortex (schematically). The vectors I (thin blue 
arrows) and M (thick red arrows; present not in scale) are depicted for the area of the core (dashed 
circle) and far from it (dotted circle); the full red dot in the origin indicates the value M = for 
the state with I perpendicular to the plane. 

where p f« 4.1 is the numerical parameter. For long cylinder with L ^> R the value of 
iV ~ 1/2 and the vortex state becomes favorable if the radius R exceeds some critical value 

Rcvit i 

'dip » I A 



R > iU - ^— j , / dip - y (8) 

where /dip determines the spatial scale corresponding to the magnetic dipole interaction. 
Note that this quantity comprises a large parameter if e /if D ~ 30 — 100, and /dip lo- 
in the case of a thin disk, L <^ R the demagnetization field energy could be revealed as 
£(homo g ) = 2nRL 2 Ml M ln(4:R/L) [9], and the vortex state is energetically favorable for 
RL > (RL) clit = 2/^p. 

For concrete estimations we take the parameters of iron borate, A =0.7T0~ 6 erg/cm, 
K =4.9-10 6 erg/cm 3 and AttM dm =120 Oe. Then we obtain that l = 3.8 nm, i.e. the 
core size is of the same order of magnitude as for typical FM (for permalloy Iq = 4.8 nm). 
The value Z^p is essentially higher, for iron borate /dip = 220 nm. Combining these data 
one finds for the long cylinder i? crit =0.9 fxm. For a thin disk sample the characteristic 
scale has submicron value, (RL) crit = 0.4 /im. The similar estimations are obtained 
for orthoferrites, and somewhat higher values for hematite. Thus, despite the fact that 
characteristic values for the dipole length /di P for FM and AFM differ hundredfold, the 
characteristic critical sizes differ not so drastically (for permalloy R CT it ~ 100-200 nm). It 
is caused by the aforementioned fact that the magnetic field created by the vortex core 
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is completely absent for the AFM vortex. The situation here is common to that for FM 
nanorings, where the vortex core is absent and the vortex state is more stable that for FM 
dots. 

Despite that the vortex core size in FM dots is rather small, the core contribution to W m 
for ferromagnetic particles of rather big radius R > 0.5 fim is negligible, but it becomes 
essential for small particles with R close to the critical size. Note as well that the vortex 
core magnetic field in the FM destroys a purely 2D distribution of M like fll]), and the core 
size changes over the thickness of the particle. For the AFM vortex the value of H m equals 
exactly zero, and truly 2D distribution of I and M, independent on a coordinate z along the 
body axis, see ([SI SI), is possible. 

Since magnon spectra of bulk FM and AFM differ significantly, one can expect an essential 
difference for magnon modes for vortex state AFM and FM particles. Remind briefly the 
properties of normal modes for disk shaped vortex state FM particles. For such samples, the 
presence of discrete spectrum of magnon modes, characterized by the principal number (the 
number of nodes) n and the azimuthal number m, is well established fla. This spectrum 
includes a single low-frequency mode of precessional motion of a vortex core (n = 0, m = 1) 
which has the frequency in subGHz region 

m = ft and also a systenn of .lightly splinted d„nb,et S with the azintnthal n um be rS 
m = ±|m|, with frequencies oo\ m \,n ^ W-\m\,ni but (jJwn\ >n ~ V-\m\,n ^ ^\m\,n |2|- The same 
classification is valid for vortices for local EP FM IllJl . Wysin had demonstrated the direct 

n 

correspondence of gyroscopic character of vortex dynamics and doublet splinting [17]. 

For an AFM vortex state particle each of two magnon branches, QFM and QAF, produce 
a set of discrete modes with given n and m, however their properties are different compared 
to that for a FM dot. The aforementioned formal Lorentz-invariance of spin dynamics of 
AFM manifests itself for motion of a AFM vortex core: the dynamical equation for the 
core coordinate X possess an inertial term, M v d 2 X /dt 2 , where the effective vortex mass 
M v = E v /c 2 17J, 118J. For this reason, the vortex core dynamics is not a precession, as 
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211 ] . but rectilinear oscillations, 



for the gyroscopic Thiele equation for FM vortices 
X{t) = a cos(u; v t+(/>o) degenerated with respect to the direction a and 4>o, with the frequency 
lu v = a/ k/M v , k determine the restoring force F = —kX for the vortex. For EP AFM model 



with W m = such dynamics has been observed by direct numerical simulations 
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22|. For 



the vortex state particle with R > R CT a the value of k is determined by the demagnetizing 
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field, k = 10 • 4ttM£ m L 2 /9R [IS], and 



2cM DM VT0T 

= — . (9) 

3y/ARln(pR/l ) 

A simple estimate gives that u v , as for FM vortex, is in subGHz region, but with different 
(approximately square root, instead of linear for FM vortex) dependence on the aspect ratio 
L/R. The other modes from this set far from the vortex core are characterized approximately 
by in-plane oscillations of / and M. As their frequencies are small, uo <C '7-Hdm> for these 
modes the magnetization M is determined mainly by the in-plane static contribution (JT]), 
and the formulae for the demagnetization field energy for FM vortices can be used. For 
these modes the frequencies are of the order of a few GHz, with approximately square root 
dependence on the aspect ratio L/R, the details will be present elsewhere. The absence of 
gyroscopical properties for the cr-ME is also manifested in the fact that for a AFM vortex 
the modes with the azimuthal numbers m = \m\ m = —\m\ are deg enerated, i.e. splitting 



of doublets with m = ±|m|, typical for the FM vortex, is absent [221 . [23] . 

For a vortex state AFM particle the high-frequency QAF branch of magnons begets a set 
of discrete modes with frequencies of the order of u> g , i.e. hundreds GHz. For these modes 
far from the vortex core oscillations of the vector / are out of plane. For their description 
the dipole interaction is not essential and the results obtained earlier for the vortex in EP 
AFM [22j, [23j can be used. The mode frequencies o; n>m are close to u g , and the difference 
w„ im — u g decreasing as the dot radius increase as c 2 /u g R 2 , with one exception: within 
the set of radially-symmetrical modes with m = a truly local mode is present, with the 
ocalization area of the order of 5/q and with the frequency u\ ~ 0.95co> g independent on R 
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The usage of QAF modes for vortex state AFM particles, particularly the truly local 
mode, would allow application of the idea of magnonincs for higher frequencies till 0.3 THz. 
A developed theory would be applied for other AFM systems like a FM bilayer dot containing 
two thin FM films with AFM interaction between them, described by the field if ex . If H ex 
is large enough, H cx > AttM s , anti-phase oscillations of magnetic moments of the layers 



produce high frequency modes with frequencies of order of J , yH ex uj™, where uj™ are the 
frequencies of modes for a single layer dot. 

To conclude, for submicron particles of typical canted AFM the ground state comprises 
topologically non-trivial spin distribution. The magnetization of each sublattices Mi and 
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M 2 are characterized by a vortex state with a standard out of plane structure, but the net 
magnetization M = Mi + M 2 form the planar vortex, where the magnetization in the vortex 
center turns to zero. Vortex state AFM particles possess a rich variety of normal magnon 
modes, from rectilinear oscillations of the vortex core position with sub-GHz frequency till 
out of plane modes with frequencies of order of hundreds GHz, including truly local mode. 
This work was supported by the grant #219-09 from Ukrainian Academy of Science. 
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